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Pressure distributions are calculated for some symmetrical wing-body 
combinations at zero lift. The theory of the calculations is based on 
the assumption of extremely slender wings and bodies and yields results 
for both subsonic and supersonic speeds. The examples considered are 
swept wings of constant chord mounted on bodies of nearly cylindrical 
form. 


Of particular interest is the effect of indenting the body on the 
distribution of pressure over the wing. When the indentation is such as 
to maintain a constant total area of the cross sections normal to the 
stream, the theoretical pressure disturbances remain small throughout the 
transonic range. With such indentation the isobars tend to remain smooth 
and nearly parallel to the sweep of the wing surface. 


INTRODUCTION 


In several papers, important extensions to the Munk -Jones slender - 
body theory (refs. 1 and 2) for lifting wings and bodies have been made 
to include the theoretical effects of thickness on the aerodynamics of 
wings and wing-body combinations. Ward (ref. 3), solving the linearized 
differential equation for the perturbation velocity potential by opera- 
tional methods, and employing asymptotic expansion of the solution, 
investigated the flow around bodies of general cross section at supersonic 
flight speeds. By a different procedure, similar results for a wing, 
body, or wing-body combination at subsonic speeds- have been developed 
by Heaslet and Lomax (refs. 4 and 5)* An analysis for subsonic flow was 
also carried out independently by Adams and Sears (ref. 6) who, in addi- 
tion, made an extension for not -so -slender wings. Confining themselves 
to wings at zero angle of attack, Keune (ref. 7) and Oswatitsch and Keune 
(ref. 8) have resently obtained a slender-body theory that is slightly 
different from those of references 4, 5, and 6. In a more recent report. 
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Harder and Klunker (ref. 9 ) have applied the basic ideas of the slender- 
body approximation to the nonlinear transonic equation for the velocity 
potential. 

The principal object of the present investigation is to apply the 
general method of reference 4 in calculating the pressure distribution 
for some special cases of nonlifting slender wing-body combinations in 
subsonic and supersonic flow. The wing of the combinations is swept back 
and has a symmetrical section with rounded leading edges. Determination 
is made of the pressure for the wing alone and for cases when the wing is 
mounted on a circular cylinder or combined with a body indented such that 
the axial variation of cross-sectional area of the combination is constant. 
The effects of Mach number and sweep angle are included in the results 
presented. 

• LIST OF IMPORTANT SYMBOLS 

bo value of x at which Sq = s(x) 

Cq root chord 

-2u 

c^ pressure coefficient^ -r;— 

^ Uo 

c*p pressure coefficient on indented wing-body combination 

t' 

Iq over -all length of the wing 

m slope of wing leading edge (See sketch (f).) 

M 
r 

Re 

Ro radius of cylindrical portion of body 

R(x) radius of indented body of revolution . 

sgn(x-g) sign of (x-|) 

So maximum value of s(x) 

s(x) local semispan 

S(x) local cross-sectional area of wing alone 
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free -stream Mach number 
polar distance in y, z plane ' 
real part of a complex quantity 
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t(x) 

spanwise distance from n 

: axis 

to wing trailing edge 

u 

perturbation velocity in 

the 

X 

direction 

Uo 

velocity of free stream 




V 

perturbation velocity in 

the 

y 

direction 

vr 

radial component of velocity 

in 

yz plane 

w 

perturbation velocity in 

the 

z 

direction 


Cartesian coordinates (x 
upward) 

downstream, y to starboard 

3 

Jw - 1 | 




e 

polar angle in yz plane 




A(x,y) 

slope of wing surface in 

X 

direction 


<p perturbation velocity potential 

<P 2 part of potential satisfying cpyy + cpgz = 0 

^ complex variable (y + iz) 

comp'lex variable (yi + izi) 

T constant related to Tq 

Tq maximum thickness of wing section 

Subscripts 

B body 

Z lower surface of the wing ( z = 0 plane) 

sub subsonic 

sup supersonic 

u upper siarface of the wing (z = 0 plane) 

W wing 
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ANALYSIS 


The thickness distribution of a symmetrical wing is prescribed. 

Under the restriction that the thickness is small and that the configura- 
tion is slender, formulas to be applied later (see "Applications" section) 
will now be briefly presented for determining the pressure coefficients 
for particular nonlifting wing-body combinations. Equations are first 
given for the wing without a body, and are then modified for cases when 
the wing is mounted on a circular cylinder or on an indented body of 
revolution. 


Nonlifting Wing Alone 


The differential equation and boundary conditions .- Expressed in 
terms of the perturbation potential cp(x,y, z), the basic linearized par- 
tial differential equation for subsonic as well as supersonic flow is the 
familiar Prandtl-Glauert equation 


(1 - + 9yy + CPzz = 0 


( 1 ) 


where Mq is the free-stream Mach number. If the surface of a wing 
z(x,y) is given, solutions to the differential equation must satisfy the 
boundary condition that the flow is parallel to the wing surface. When 
the wing is thin, it is sufficient to satisfy this requirement in the 
plane z = 0. Analytically, the expression of the condition is 


Sz 


= wu(x,y) = Uo = UoXu(x,y) 


z=+0 


( 2 ) 


where Uq is the free-stream velocity and w^(x,y) is the vertical 
Induced velocity on the upper side of the z = 0 plane. 

Velocity potential .- When the flow is supersonic, the formula for 
the perturbation potential subject to the boundary condition (2) is known 
to be (e.g., see ref. 10) 


cp(x,y,z) = ~ n J'J' arccosh 

T 


py (y - Tl)‘ 


dll (3) 


+ z‘ 
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where the region of integration T is the portion of the plan form lying 
within the Mach forecone from the point In the case of subsonic 

flow fields, the solution of the differential equation (l) may be expressed 
in the form 


9(x,y,z) 





arcsinh 


— — =^=d| dTl 

3.y(y - r\)^ + 




with T now extending over the entire plan form. These two solutions 
yield the potential due to a distribution of sources of strength propor- 
tional to the slope of the wing surface ^u(^>y) • 


When the wing is slender, that is, if 


JL 




is considered 


X - I 

very small, further approximations to the linearized potential (3) and (4) 
can be readily made. In this event, one may employ the approximate 
relations 


arccosh 


X - 5 


3 J{y - + 


^ In 


2(x - I) 

iy(y - + 


> (5) 


arcsinh 


X - 5 2 |x - ^ 

• ^ sgn(x - 5 ) In 


3^(y - + 


3-s/(y - T))^ + 




where the symbol 
taken. 


sgn(x - i) means that the sign of (x - 5) is to be 


Consider now a thin pointed wing 
of symmetrical section with straight 
or sweptforward trailing edges as is 
shown in sketch (a) , Use of the 
relations (5) in equations (3) and (h) 
then gives the result (ref. 4) 


cp(x,y,z) = cp2(x,y,z) + g(x) 


( 6 ) 
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where 


cp2(x,y,z) 


^ + z^ldT] 

S2(x) 


( 7 ) 


and 


g(x) 


H2 A 

2« dx 



2(x - 0 
3 


d| 


(8a) 


for supersonic flow, and 

X Un 8 r^° 2 lx - I 1 

' Tw ®Sn(x - 4)S»(|)ln ^ d5 (8b) 

o 

for subsonic flow. The function S'(x), the derivative of the cross- 
sectional area of the wing in a yz plane, is found from 


S'(x) 



s^{x) 


( 9 ) 


It is seen from equation (6) that the slender -body approximation to 
the linearized supersonic and subsonic potentials consists of two parts. 
The first part (iPg) in each case is independent of Mach member and is a 
harmonic function in the transverse plane; that is, it satisfies the two- 
dimensional Laplace equation 




( 10 ) 


The second part (g) depends on the cross-sectional area of the wing and 
is a function of x and P only. 

Inspection of equations (6) and (8) shows that the value of the 
potential for a particular wing at some Mach number Mi can be wi'itten 
in terms of the potential given at another Mach number Mq. Thus 


CONFIDENTIAL 



NACA RM A5itJ07 


CONFIDENTIAL 


7 


*?sur)n Oir S*(x)ln 




Uo 3i 

'Psubp^ = ^suh^ 2rt ^ 


(lla) 


with Pi = y|Mi^ - l| and 3o = “ 1 1 • The supersonic and subsonic 

potentials may also be related. If one has already found an expression 
for the supersonic potential, he can then obtain the subsonic potential 
by means of the equation 


U. 


9, 


sub 


= ”sup ST S'(x)ln 


Psuh 

P 


sup 


Uo r 


s» (I) 
I 


X 


d| 


(IIL) 


simplified theory, provided the 
\Uo 


Certain symmetrical wings whose trailing edges are 
sketch (h) can also be treated by the 
chordwise variation in shape of the 
cross sections is sufficiently smooth 
and gradual that the assumption of 
two-dimensional transverse flow may 
reasonably be applied. Since the 
wing is at zero angle of attack and 
the flow is symmetrical about the 
xy plane, the source distribution 
^u(x,y) in the plane of the wing in 
region W between the axis and the 
trailing edge is set equal to zero. 

If the limits of integration are 
properly adjusted, equation (6) will 
then formally still apply and yield 
expressions for the potential in the 
two regions 1 and 2 of the sketch. 

Thus for region 1, in the z = 0 plan 
one obtains 


swept back as in 



“^sup 


jt 


s 

Au(x,il)lnjy^ - T]^|dTi " 


o 


Uo 

2^ Sx 


Si«(|)ln I (x - ^)d|, 
o 


(x < Cq) 


(12a) 
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and in region 2 


Uq rs Uq 3 -Ro 

9sup = — J Xu(x,Ti)in |y^ - Si»(|)in(x - |)d| - 




‘’‘-'o’ 


(12b) 


where 


s(x) 


Si»(x) - ^ z^(x,n)<iTi, 0 < X < Co 


( 13 ) 


s(x) 


S 2 *(x) = ^ ^ / Zu(x,Ti)dn, 


Co < X < I, 


t(x) 


The subsonic solution is obtained from the supersonic solution by using 
equation (lib). 

Pressure coefficient .- After determining the potential for a wing 
from the equations in the foregoing section, the pressure on the surface 
of the wing is found by differentiation. The pressure coefficient is 
related to the perturbation velocities by the equation (ref. 4) 


Cp(x,y,0) 


_Uoax‘"u7vW ^ Uo^ 



(I4a) 


which is an approximation to the complete Bernoulli equation consistent 
with both the linearized differential equation (l) and the assumption of 
slenderness. For planar problems, further simplification achieved by 
neglecting the nonlinear terms (S<p/Sy)^. and (Scp/Sz)^ yields satisfactory 
estimates for the coefficient. 
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Thus 


Cp(x,y,0) 


1 . 


(14b) 


where <P is obtained from equations (lib) and (12). 

It should be pointed out, however, that equations (lib), (12), and 
(l4) will not furnish realistic results for the pressures^ for a wing in 
either subsonic or supersonic flow unless certain restrictions are imposed 
on the gradient of cross-sectional area and its derivative S"(x) . Con- 
sider, for example, the pressure found along the line AB where regions 
1 and 2 in sketch (b) join. In supersonic flow, one finds from equa- 
tions (12) and (l4b) that at this line there is a jump Acp in the pres- 
sure coefficient given by 


= -P, - -Pa = i e^O 

Pv s2(co) - y2 


S ’'u(x.t) 


Cr 


in 

+ e PyJ 


(15) 


where 


AS"(co) = Si"(co) - S2"(co) 


and where the usual assumption that S'(0) = S*(cq) = 0 has been made. 
Evidently, along the line AB a logarithmic singularity will occur in 
the pressure distribution in going from region 1 to region 2 if 
Si"(co) / S 2 ”(co).' (The infinity would of course be higher than loga- 
rithmic if AS"(co) 1® singular.) The formulas for the pressure coef- 
ficients presented in this report and in reference 4 are therefore good 
only ^or cases where the plan form and the slope of the wing surface are 
sufficiently smooth so that there are no abrupt changes in either S* or 
S". (This holds true also for the formulas given in references 3, 6, and 
7 . The restriction on S,", however, may be somewhat relaxed in employing 
the slender-body theory of reference 8.) 

^Although the theory may give spurious infinite pressures on certain 
portions of the wing, the results obtained for the wave drag by integrat- 
ing the product of pressure and surface slope over the wing may be finite 
and reasonable. 
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Even if there is no jump in the pressure coefficient at the line AB, 
the slope of the pressure curve may not be continuous. It can be shown, 
for instance, that a singularity will occur in the slope of the curve for 
the pressure coefficient in case S"'(co) is discontinuous. 


Wing on a Cylindrical Body 


The 
will now 


equations presented in the foregoing sections for the wing alone 
be modified to yield formulas for calculating the pressure dis- 
tribution on a combination composed of a 
symmetrical wing zu(x,y) = -zj(x,y) 
mounted on an infinite circular cylinder 
having radius Ro(see sketch (c)). The 
surface of the wing chosen here will also 
be considered symmetrical about the xz 
plane, but the method applies equally well 
if this is not the case. The procedure 
followed is essentially the same if the 
fuselage is any body of revolution instead 
of a circular cylinder. 




C, -plane 


!Z, 


■s. -R, R, s, 


Sketch (c) 




In studying such combinations, it is 
usually convenient to introduce a second 
coordinate system. Let the yz plane be 
represented by a complex variable 


.16 


C = y + iz = re" 
and then consider a plane 

= yi + izi = rie^^i 


obtained from the C plane by the 
Joukowski transformaWon 


(16) 




( 17 ) 


The transformation (17) maps the C plane onto the plane so 

that a circle representing a section of the body in the C plane is 
mapped onto a portion of the real axis in the plane, while the part 

of the real axis outside the circle is transformed into an adjoining part 
of the real axis of the plane. (See sketch (c).) It can easily be 
shown from the equation that the geometric relations 
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Ri = 2Rc 


yi = y + 


y ^ 


yi^ > Ri' 


yi = 2RoCOs e, yi^ S Ri‘ 


(18) 


Ro" 

Si = s + -f- 


hold for zi equal to zero. 

The technique employed is to transform the boundary conditions from 
the ^ plane to the plane and to find a solution in the latter 

plane. The solution is then transformed back to the physical plane for 
the completion of the problem. 


Perturbation velocity potential at the wing surface .- A consequence 
of the conformal transformation is that the complex velocities in the 
two planes ^ and are related by the equation 


V - iw = (vi - iwi) 




(19) 


or, in polar coordinates (ref. 11), 


V = Vi 


1 - 


Rr 


COS 29 






+ j sin 20 


w = Wi 1 - 


— ) cos 20 


- vi 


R, 


sin. 20 


> (19a) 




Vp = vi cos 0 + Wi sin 0 


"R 1 

1- 1 — 




+ 2wi I -p-j sin e 


\ 


J 


From these equations it follows that the boundary conditions in the zi = 0 
plane are 
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Vr ^ Uo(dRo/dx) 
2 sin 9 2 sin 9 


yi^ < Ri^ 


> ( 20 ) 


Wi 


V 

1 - Ro^/y^ 


Ri^ < yi^ < si^ 


with w related to the streamwise slope of the wing surface by equa- 
tion (2) . 


The two-dimensional solution Cp 2 in the plane at zi = 0 is, 

from equation (J), 




Putting now 


<P2(x^yi^0) = wi(x,Tii)ln|yi - TiildTii 

-Si 


( 21 ) 


Til = q + 


Ro^ 


and using relations (l8) and the boundary conditions (20), one . finally 
has for equation (21) 


T 1 

Uo r 

q>2(x^y) =~ A^(x,Ti)ln 


Rr 


(y^ - Ti^)(y^n^ - Rq^) 

y^T,^ 


dT] 


(22) 


If one of the functions g given by formulas (8) is then added to this 
equation for cpg , the velocity potential on the wing may be obtained for 
either supersonic or subsonic flow. Thus, for the supersonic perturbation 
potential on the upper surface, there results 




sup 


s(x) 

U_ r I 

= — / 


(y^ - Ti^)(y^Ti2 - rq4) 


2 2 
y n 


dTl - 


^ A /" S«(l) Inf (x 


I)d4 


( 23 ) 


and for the subsonic potential 
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*Psub 



(y^ - - Ro^) 


dTi - 


H° i- 

ij-Jt Sx 


r 

J sgn(x - |)S»(|)ln I I X - I |d| 


i2k) 


Pressure coefficient at the wing surface #- The pressure coefficients 
on the surface of the wing of the combination are obtained by differenti- 
ating equations ( 23 ) and (24) with respect to x and then employing 
relation^ (l4b) . Thus 


^sup 


s(x) 


- Ro‘) 


Ro 


2 2 

y Ti 


dTJ + 


«&/ S'«)lnf(x-t)dl 

o 


( 25 a) 


and 


"Psub 


s(x) 

I-/ Au(*.l)ln 

Ro. 


(y^ - Ti^)(y^n^ - Ro*^) 

7? 


di] + 


1 


2^^ / sgn(x - |)S»U)ln ^ |x -|| d| 


(25b) 


or, with the aid of equation (lib). 


"Psub 


L 

1 ^sub 1 ^ r 
- ^Psup - « S (x)ln J 


S'(g) 

X - i 


d| 


( 26 ) 


2 

As in the case of planar problems, the squared terms in the pres- 
sure relation can be neglected in considering a combination whose body 
is a circular cylinder. (See ref. 4.) 
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The coefficient for a particular wing at some Mach number Mi, written 
in terms of the coefficient at another Mach number Mq, is 


supp 


= c. 


sup 


Po 


+ i S"(x)l. ^ 


(27a) 


and 


sub 


'Pi 


and 

Po = 




(27b) 


It should be noted that formulas (25a) and (25b) are not uniformly 


valid for all values of Er 


\Uo 



If the radius Rq is equal to zero, the 
• equations reduce immediately to those 
given previously for the wing without 
a body. As the radius approaches 
infinity, however, one finds that for- 
mulas (25a) and (25b) yield a value of 
the pressure coefficient that differ 
from the results for the wing alone by 
S”fx) 1 

a term equal to — ^ In ^ • Employing 

the equations for values of Rq that 
are large violates the assumption of 
slenderness according to the approximate 
theory used herein. 

When the wing mounted on the com- 
bination is swept back, as shown in 
sketch (d), the slope Au(x,t]) is taken 
equal to zero in the gap between the 
trailing edge and the body. The for- 
mulas obtained in regions 1, 2, and 3 
for the supersonic pressure coefficient 
on the upper surface of the wing are 
then as follows: 


Region 1: 


2 a 


Hr 


(y^ 


i^)(yV 


y^Ti‘ 


Ro^) 


dTi + 


1 

n ^ / Si»(Oln f (x - 5 )de, (0 < x < co) 


(28a) 
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Region 2: 

s(x) 

=Pa = ■ f I;/ ^u(>:.n)ln 

t(x) 


(y^ - - Rq^) 

y^T)® 


drj + 




~ k ipf S 2 »(|)ln I (x - 

o eo 

(co < X < bo) (28b) 


Region 3: 


So 


'P3 = - I 8^ / ^u(^.Ti)ln 

t(x) 


(y^ - n^)(y^n^ - Ro^) 


2 2 
y Ti^ 


dT] + 


«o °o 

Si'(e)ln(x- - |)d? + i y S 2 ‘(|)ln(x - |)d^ + 

Co 


X 

I § / Sa’COln I (x - Ode, (bo < x < l^) 


( 28 c) 


where bo is the value of x at which Sq = s(x) . Use of the relation 
( 26 ) in conjunction with the above three formulas will furnish equations 
for the pressure coefficient in subsonic flow in the various regions. 


Wing on an Indented Body 


In the previous formulas for the pressure distribution attributable 
to thickness, the coefficient becomes infinite^ when the Mach number 
approaches unity because the formulas contain a term involving S*(x)ln — . 

It is thus possible to construct slender configurations which will give 
a theoretically finite pressure coefficient even at the speed of sound. 
Under the assumptions of slender-body theory, this may be achieved if the 
gradient of cross-sectional area of the configuration in a yz plane 
vanishes identically. Combinations constructed in this manner on an 

This of course is not a property of slender-body theory alone. 

Except in particular cases, steady-state linearized theory also in general 
yields infinite pressures at g = 0. 


CONFIDENTIAL 



16 


CONFIDENTIAL 


NACA RM A 5 UJO 7 


infinitely long fuselage are slightly indented in the vicinity of the 
wing and possess, according to the theory, the important property of 
having zero wave drag. 


Velocity potential at the wing surface .- Let the surface 
Zu(x,y) = -zj(x,y) of a symmetrical wing be specified and assume that 
the indented fuselage is a body of revolution which deviates slightly 
from a basic circular cylinder of radius Rq . The cross-sectional area 
S(x) of the exposed wing is then given by 

po 



S(x) = h ■ z^^(x,Ti)dT] 
R(x) 


(29) 


where R is the radius of the indented 
body. (See sketch (e).) In order for 
the streamwise gradient of cross- 
sectional area of the entire combina- 
tion to be zero, the relation 


flR^(x) = «Rq 2 - S(x), S(x)<<flRo^ 


. -S'(x) 


(30a) 


(30b) 


must hold. Since the quantity S(x) is known from formula ( 29 ), equa- 
tion (30a) can thus be solved for the radius R of the body. 

Now the perturbation potential in the plane at zi = 0 has the 

form 


1 

=- Wi(x,Tii)ln I yi - rji j diii 
-Si 


(31a) 


or 


1 1 r"' 

=-/ Wi(x,Tii)ln|yi - Tii|dTH + ^/ wi(x,t 1 i) Injyi^ - Tii^ldiu 

(31b) 
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with 


Si = s + — , 


r2 

7i = 7 + — , = 2R 


(32) 


The term corresponding to the g(x) of equations (8), which is to be added 
to this cp2^ is zero since there is no change in the cross-sectional area 
of the specified shape of the yz plane. 


From the last relation in equation (l9a) it follows (upon putting R 
in place of Rq) that in the interval < ri^ tne boxindary condition 

.may be expressed 


vr _ (dR/dx)Uo _ 

2 sin 0 ■ 2v/l - (Tii/ri)2 ’ 


(33) 


For very thin wings, and for deviation Rq - R of the same order as the 
wing thickness, the condition in the interval ri^ < T|j ^2 ^ approxi- 

mately 


V _ Uq(5zu/^x) 

1 - E^/yf " 1 - r 2 / ti 2 


(34) 


Making use of equation (33)> one can then write the first term on 
the right in equation (31b) as 


ri 

Wi(x,T]j^)ln ly^^ - n^^ldq^ = UqR ^ In y, y^ > R^ (35) 

-^1 

where^ from equation (30b), the quantity R(dR/dx) may be replaced by 
-l/2jt times the gradient of cross sectional area S*(x) of the wing. 

In the second integral in equation (31b), set 


ni = 



and use relations (3^)* The final form for the velocity potential at the 
upper surface of the wing thus becomes 

4 >(x,y,0) = - ^ S'(x)ln y + 

I 


Uo 

— / Au(x,q)ln 


(y^ - T]^)(y^Ti^ - R^) 


y^T)2 




y > R 


(36) 
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Formula ( 36 ) for the velocity potential for a slender combination 
whose body is indented according to the area rule (ref. 12) is evidently 
independent of Mach number and therefore holds (within the assumptions 
of the simplified theory) for subsonic, transonic, and supersonic speeds. 
If the limits of integration are adjusted, the equation may also be used 
for combinations with sweptback wings . 

Pressure coefficient at the wing surface .- On the surface of the wing 
the pressirre coefficient can be found from relation (l4a), that is. 


c 


Pw 


"2_ 8cp 1 


_Uo 8x Uq2 

\8yy 



(37) 


where the potential 


cp is obtained from equation ( 36 ). 


APPLICATIONS 


The formulas which were given the preceeding part of this report will 
now be applied for the purpose of performing detailed calculations of the 
pressure distribution for some particular nonlifting combinations having 
a symmetrical wing with constant chord. 


Wing Mounted on a Body of Circular Cross Section 


Pressure coefficient on the surface of the wing .- The wing of the 
combination"^ considered here is mounted on a circular cylinder and has 
an upper surface defined by 

Zu(x,y) = T(com - s + y)N/{s~^^~yy(c^m~rVTllO (38) 


with 


s = mx + Rq 

zu(x^y) = Zu(^^-y) = -zi(x,y) = -Z2(x,-y) 


and 


8tq 

T = 

3m2coW3~ 


( 39 ) 


■*This combination will be referred to as the basic combination. 
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the constant Tq being the maximum thickness . The profile resembles a 
Joukowski section and is the same for all values of y. On the starboard 
side of the wing, the slope of the surface is given by the relation 


Sx 


Au(x,y) 


Wu(x,y) 

Uo 


mr 

2 


(c^m - ks + 


4y) 



s + y 

y 


m 


(klb) 


From this equation it is seen that the slope possesses a square-root 
infinity at the leading edge and is zero at the trailing edge* 

The gradient of cross-sectional area of the wing is 
Si*(x) = 4Tm^(c^ - x)'n/x(cq - x) = 4m*z^(x^Ro)> (O < x < co) (^la) 

in region 1 in sketch (f) and is 

S 2 '(x) =0, Co < X < (so - Ro)/m (klb) 

in region 2; for region 3 > 

83 ’ (x) = -4m Zu(x,Sq), (sq - Rq)/^ < X < (sq - Rq + Com)/m 

Since S" is not continuous across | ^ 

the line AB, the approximate theory 
eu^loyed here will give unrealistic 
results in the regions of the tips . 

In calculating the pressures for the 
wing-body combination, attention will 
therefore be confined to regions 1 
and 2 only. The calculations to be 
given in subsonic flow will be based 
on the assumption that, the tips are 
located far downstream and have no 
effect on the other two regions. 

Use of the equations (4o) and 
(4l) in formulas ( 28 a) and ( 28 b) yields 
the following final results for the 
supersonic pressure coefficient on the 
upper surface of the wing for regions 
1 and 2: 
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Region 1: 


"Psup 


= l6R„tan~^ / ' + 2m(co - ^x) / 

rt|_o Vcq-x V 


Co - X 8xy(co - x) 


In 


3co(y^ - Rq^) 


F(x,y) + F(x,-y) - 2 F(x, 0) + F^x,-^+ f(x, - 

(O < X < Cq) 


ih2a) 


Region 2: 




'Psup 


/ X - Co , /s - Com 

8Ro + 2ni(4x - Cq) J — + 2(4s - Cotny^ ; + 


y s “ C ! 

— TT 


s - CqUI + y / Rq \ /s - Com - Ro^/y 

y V y y*/ s - Ro^/y 


c^m - 4s - 4 y 


Ro^\ / s - Com + Ro^/y 


s + Ro^/y 

(cq < X, s - Com <y<x; s = mx + Rq) 

(42b) 


where 

/ ^ / > IS / s - Com + 

F(x,t)) = (com - 4s - kr\) 

V S + T1 


[f- 


tan 


_i (co - x)(2s + Ti) _ (tix + RoCo) 


2 n/x(Cq - x)(s + T])(s - Com + Ti) 


Cq) 1 

Tn)J 


(43a) 


when c^tn - s - t] < 0, and 
F(x,ti) = (com - 4s - 4T1) 


CqTII - S - T] 


S + T] 


when 


In 


Com - 


(co - x)(2s + Ti) _ (tix + RoCo) + ^/x(co - x)(s + 'n)(com - s - t]) 

(43b) 

s - n > 0 
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It is seen from equation (26) that the formulas for the pressure 
coefficient in subsonic flow can be obtained by adding two terms to the 
above equations for the supersonic pressure coefficient. Since the value 
of S' in region 2 is identically zero, the terms to be added for region 
1 become 

Co 

■ ^ Jc y' X -^1 " Tm^[4x - 3cq], (x < Co) (44a) 

o 


and for region 2 


2.n 6x 



d| 


‘ k 


rm^ 4x - 3 cq + ( 


Co - 



(x > Co) 

(44b) 


The subsonic pressure coefficient on the upper surface of the wing can 
therefore be written in region 1 as 


•sub 


"Psup 


^ fl(4x - 3 Cq) + 2(4x - Cq)J ~ ^ In (45a) 


and in region 2 as 


■ sub 


c + m^T 4x - 3 Cq 

^sup L 


+ (c 


o “ 



(45b) 


where c„ 

*^sup 


is given by equations (42) , 


Formulas (42) and (45), which give the pressure coefficients on the 
wing of the combination indicated in sketch (f), will be plotted and the 
results discussed in a later section. An obvious result noted now is 
that in region 2 the formulas do not depend on Mach number for either 
supersonic or subsonic flow, but that the subsonic coefficient in this 
region is always greater than the supersonic coefficient because the 
inequality 


m^T 


4x - 3 cq + (cc 



( 46 ) 


holds for all values x ^ Cq- 
1 
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As a check, let us nov consider the asyn5)totic behavior of equa- 
tions (38), (k2b), and (45b) far outboard along the wing. For this pur- 
pose, we introduce a change in coordinate system by means of the trans- 
formation 


X = xnsin 0 + ynsin 0 

y - Rq = -xncos 0 + ynSin 0 


(47a) 


\Uo 


and 



(See sketch (g).) Setting tan 0 = m, 
one can then write 


X = 


y - Ro = 


iPXn + yn 
■Jl + m2 

myn - Xn 


yn = 


n/i + m2 
mx - (y-Ro) 


>(4Vb) 


n/1 + 

X + m(y- Rq) 
■Jl + m^ 




Cj^ = c^m/s/l + m^ 


= m/Tl+m? 

Un = UQm/*yi + m^ 




Pn = l/7l + m2 


(48) 


Equation (38) thus becomes 


^ — (cn - 
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which agrees with the approximate equation employed in two-dimensional 
section theory for a thin Joiikowski base profile of thickness ratio 

— — • (See ref. 13 •) The pressure coefficient for such a two- 

3*/3c„ 

dimensional airfoil is 


2 dcp 


8^0 ■ 
3'/3cn^n - 


3 - 



or, written in terms of the xy system. 


(^9a) 




8tq 

3‘/3co^ 


[3cpm - lf(mx - y + Rq) ] 


(49b) 


This equation is in agreement with the asymptotic expression of the pres- 
sure coefficient for region 2 obtained far outboard along the wing from 
either equation (42b) or (45b). 

Pressure coefficient on the body .- The formula for determining the 
pressure coefficient on the surface of the body is (ref. 4) 




(50a) 


When the body is a circular cylinder, this relation reduces to 



(50b) 


It is apparent that application of equation (50b) requires a knowledge 
of the value of the potential q) in space, so that formulas ( 28 ), which 
hold for z = 0, must be "analytically continued." One way that will in 
effect accomplish this is as follows: 

Set 


.. . „ ie 

^ = y + iz = R^e 


and form the function 

^(x,C^0) = cp(x,^,0) + ii|r(x,^,0) (51) 
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where cp(x,^,0) is taken from equations ( 23 ) and (24), with the function 
defined by^ 


s(x) 

i|f(x,^,0) =J^ Wu(x,q)dTi ( 52 ) 

The constructed function ( 5 I) will obviously satisfy Laplace's equation 
in the yz plane; and it can be verified without difficulty that the 
boundary condition Vp = 0 is also satisfied. The coefficient on the 
body can then be calculated from the formula 


c 


P 




) 


where the symbol Re means that the real part 


is to be taken. 


( 53 ) 


For region 1, the final result obtained for the pressure coefficient 
on the body in supersonic flow for the basic combination considered in 
this part of the report is 


• sup 


m?T 

n 


16 Rq tan' 


‘7 


i0N 


Co - X 


- 2F(x,0) + F(x,Roe ^ + 


F(x,Roe"^®) + F(x,-Ro6^®) 


+ F(x,-Roe"^®) + 


Re 




8x(cq - x)e 
pco(e^^® - l)Ro 


ijt(cQm - 4s + 4RQe^®) 


'cr,m - s + Rr>e 


ie 


s - Roe 


ie 


(x < Cq). 
(54a) 


where the function F(x,t]) is given by equations (43) . The result found 
for the supersonic pressure coefficient in region 2 can be put in the form 

^An arbitrary real function E(x) may, in general, be added to the 
right side of equation (52). Such a function, however, would in no way 
affect the real part of S$/5x, since iE is pvirely Imaginaxy. 
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Cp = m2 T 

^sup 




8Rq + 2m(4x - Co) I • ' y + 2(4s - Com) J — 


(com - 4s 


/ D 

i0^ /s + Ro® “ Com 
4Ro 6 ) / Trt + 


s + R^e 


i0 


(com - 4s - 4Roe"^®) / + 


S + Rq^ 


-i0 


±9 


( CqUI - 4s + 4 Rq0 ) 


±e^ r " .^oe - Corn 


s - RqC 


ie 


(com ■- 4s + 4Roe~^®) 


's - Rr,e' 


•i0 


- Com 


s - R^e~^® 


(x > Co) (5^b) 


Corresponding formulas for the pressure coefficient on the body in sub- 
sonic flow are obtained from the above two equations by adding to them 
the terms on the right of equations (44a) and (44b), respectively. For 
any fixed value of y in region 2, relations (45b) and (46) furnish the 
inequality 


^Psub 


> c 


Psup 


(55) 


which, as mentioned before, is also true for the coefficients on the wing. 


Wing Mounted on an Indented Body 


The constant -chord wing whose thickness distribution on the starboard 
side is given by 


2u(x,y) = T(com - s + y) n/ (s - y)(cQm - s + y) , R(x) < y < s ] 

S = TIIX + Rq 
8 tq 
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will now be combined with an indented body of revolution which deviates 
slightly from a basic cylinder in the manner that makes the local cross- 

sectional area of the configuration 
\Uo *^oiistant. The wing has the same 

\ profile along EE as the one in the 

previous example. (See sketch (h).) 
Since the gradient of cross-sectional 
area of the combination is set equal 
to zero, the wave drag, according to 
the slender-body theory, is also zero. 

The cross-sectional axea of the 
wing in region 2 is equal to the con- 



stant 


SgCx) =4 y Zu(x,Ti)dTi = — 




(57) 


s-CQm 


and the radius of the indented body 
in this region is therefore the con- 
stant 


R.. 


“ 


Cq^ih® 


(58) 


Corresponding equations for region 1 are 


Si(x) = ^ J ^ |3m®Co3cos"i 


Coin - 2s + 2Ri 


Cotn 


Rn 


2[2(s - Ri)(7cQm - 4s + 4 Ri) - (s - Ri)(com - s + Ri)j-(59) 


and 


Ri(x) 



O<R2<R1<R0 


(60) 
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Equation (60), however, does not give an explicit value of the radius 
in region 1 because the function itself involves Ri. Moreover, 

the function is a transcendental function of Ri, so that equa- 

tion '( 6 o) can only be solved by graphical or other methods of approxima- 
tion. By means of the mean-value theorem, we can write the equation in 
the form 

Si(x) = Sii(x) + 4zu(x,d)(Ro - Ri) ( 6 l) 


where d is a certain value in the interval Ri < d < Rq, and 


, , m3r 
Sii(x) = — 




(62) 


is the cross-sectional area of the wing in region 1 neglecting the addi- 
tional area exposed by the indentation. Since the deviation Rq - Ri is 
assumed to be of the same order as the wing thickness, neglect of the 
term 4z^(x,d)(Ro - Ri) will evidently introduce only an error of the 
second order in thickness. For very thin wings, the radius Ri may thus 
be approximated by 

which is now an explicit fimction 
of X. 


Sketch (i) indicates the varia- 
tion of the radius as a function of 
x/cq for To/mco equal to 0.1 and 
values of Ro/mco equal to 0.5 and 

1 . 0 . 


Relations (l^b) and (53) can be 
immediately employed for determining 
the pressure coefficient on the sur- 
face of the wing and on the body in 
region 2 , because the portion of the body in that region is a circular 
cylinder. Use of these relations with equation ( 36 ), (5?), and ( 58 ), then 
finally yields for the coefficient on the upper siarface of the wing 

Cp = m 2 r jacks - Sc^m) + 2 (ks - + (com -ks - ky) ' + 



Sketch (i) 


c^m - ks + k 



s - Cpm -Rg^/y 


- Ha^/y 


Cotn 



s - Cpm +Rg^/y 


+ R2^/y 


- Corn < y < s) ( 6 k) 
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and for the coefficient on the body 


=n^T 


2(4s - 3coin) +2(4 s - Com) +(0^111 - 4s - 4Rge^^) / ^ 

s - Rae^® 


(cotn - 4s - 4R2e"^®)y^— +(com - 4s + 4R26^®) J 


s - Cq® - Rge 
8 - Raei® 


15 




^±G 


(65) 


These equations are employed for 


X ^ ^o - m 


Ho . 


Cq%1T 


m2 4 


the equality sign giving the value of x at the trailing-edge fuselage 
junction. The two formulas do not involve Mach number and apply for sub- 
sonic^ transonic^ and supersonic flight speeds. Comparison of formula 
(64) with equations (42b) and (45b) shows that for a fixed y in region 2 
the inequality 





> c- 


’^sub 


> Cp 


sup 


( 66 ) 


is satisfied when x > Cq . Values calculated in region 2 for the subsonic 
and supersonic pressure coefficients along a section on the surface of the 
wing of the basic combination will thus be less than the values obtained 
for the coefficient Cp on the wing surface of the indented combination i 

In region 1, the body is not cylindrical but it is found that the 
squared terms in the pressure relation (3?) may be neglected since they 
contribute only quantities involving the second and higher order in wing 

2 dffi 

thickness that are small in comparison with - 77— . Even the first 

bo ox 

term in relation (37) gives rise, for the particular combination consid- 
ered, to some small terms of the second order. Such quantities, however, 
are also fo\md to be negligible. The formulas to be presented here for 
the pressure coefficients \ in region 1 will therefore contain (like those 
for region 2) only terms of the same order as the thickness of the wing. 
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Use of equations (l4b), ( 36 ), (53), (59), and ( 60 ) thus furnish the following final expres- 
sions for the pressure coefficients on the wing surface and on the body in region 1 of the combina- 
tion in sketch (h); 


'Pi 


m^T 

Jt 


2 (mCo - 4s + 4 r 


■>/ 


Cpm - s + Ri 
s — R 1 


In 


~ 8n/ (s -Ri)(c m - s +Ri) + 




4(4s -3Com)tan-'-y^^^I^+Fi(x,y) +Fi(x,-y) -2Fi(x,0) +F^ + (x,- 5^^ 


(67) 


and 


V- _ m^T 

'Pi = IT 


-dj{s -Ri)(com- s +Ri) +4(4s -3Com)tan”^/ — + Fi(x,Rie^®) +Fi(x,Rie"^®) - 

^ V Cpm - s + Ri -i-' ^ 1 


2Fi(x,0) +Fi(x,-Rie^®) + Fi(x, -Rie’^®) +Re 


2 (cpm- 4s + 4 Ri) / In 

V S - 


esi9 




L - s + R^e 
Riei® 


ijt(cQm - 4s 
where for Cpm - s + >0 

Fi(x,ti) = (com-4s + 4Ti) 

V s - n 


( 68 ) 


In - 


mcp(Ti - Ri) 


•2s ( Cpm - s) + (2s - CQm)(T| + Ri) - 2Ri + 2-7 (s - tj) (cpm- s + q) (s - R^) (cpm - s + R 


1). 


and for c^ra - s + q < 0 
Fi(x,t]) = (com - ks + k)J 


s - Cpm - T] 
S - T] 


tan"i 2 s(cpm- s) + (2s - Cpm)(Tl +Ri) - 2TlRi 

2^j is -Ti)(s - Cpm - ti)(s - Ri)(c m - s + Ri) 


(69) 


(70) 
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NUMERICAL RESULTS AND DISCUSSION 


Graphs of the pressure coefficients for a wing alone and for the 
two wing-body combinations in sketches (f) and (h) are shown in figures 1, 
2, and 3- Figures 1 and 2 give plots for the coefficients along several 
spanwise stations of the basic combination (Ro/cQm = O. 5 ) for subsonic 
and supersonic flow at mp = O. 5 , while figure 3 presents plots of the 
coefficient for the indented wing-body combination which are independent 
of Mach number. The plot in figure k gives the variation of the pressirre 
on the body in the plane of the wing. 

Figures 1 and 2 show that, except at the wing-fuselage juncture, the 
pressure coefficients have a finite negative® value at the leading edge 
and increase to finite positive values at the trailing edge. On sections 
which are cut by the plane x = Cq (i.e., passing through the trailing- 
edge fuselage juncture) the slopes of the curves are discontinuous . It 
is also apparent that the eff.ect of the presence of the body on the coef- 
ficient does not extend very far downstream beyond the trailing-edge 
fuselage junctiire; the difference between the calculations for the wing 
on the combination and those for the wing alone, for instance, are too 
small to show up in the plots for sections more than one chord length 
from the body. 

In figure 3^ the pressiire coefficients for nearly all sections on the 
wing of the indented combination are very close to the curve for the two- 
dimensional wing. The discontinuity in the slope of the curve that was 
quite noticeable along sections such as BB of figures 1 and 2 is far less 
apparent in figure 3* 

The graph in figure 5 is the fvmction which, according to equa- 
tions ( 27 ), can be added to the values given in figures 1 and 2 for the 
subsonic or supersonic pressure coefficient along the section AA and BB 
of the wing to yield values of the coefficient at other Mach numbers. 
(Sections CC, DD, and EE lie in region 2 where the coefficients do not 

depend on Mach number.) For example, the pressure at the 

(To/co)m 

leading edge along BB in subsonic flow for mp = 0.4 is, using figures 1 
and 5 in conjunction with equations ( 27 ), -4.1 -0.2 = -4.3- 

Figure 6 shows isobaric charts of the pressure coefficients for the 
wing without body in subso nic and supersonic flow at mp = 0.^, and also 

® The occurrence of a negative pressure at the nose of a Joukowski- 
like section is the result of the thin-airfoil simplification. A more 
accurate theory would show a small region of positive pressure (i.e., 
a stagnation point with the maximum value equal to the impact pressure 
of the component of stream velocity normal to the edge) . 
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a chart giving the two-dimensional results. The figiares 6(a) and 6(b) 
indicate that the pattern of the isobars in the region behind the trailing- 
edge fuselage juncture is essentially the same as the two-dimensional in 
figure 6(c), but that in the region near and upstream of the jimcture, a 
marked deviation from the straight isobars of the two-dimensional case is 
evident . 

Isobaric maps are also shown in figure 7 for the basic wing-body 
combination in subsonic and supersonic flow, and for the indented wing- 
body combination. Figures 7(a) and 7(b) illustrate that in the region 
downstream from the trailing-edge fuselage junct\ire (even near the juncture 
itself) the isobars are not much different from those for the wing alone, 
and that the body therefore has little effect in this region. In the 
region adjacent to the body, the pattern of the isobars is qualitatively 
similar to the case for the wing alone but the pressures are lower. A 
remarkable difference between the chart in figure 7(c) for the indented 
combination and those in figiires 7(a) and 7(b) for the basic combination 
is that the isobaric pattern on the wing for the indented combination is 
essentially two-dimensional over practically the whole wing. 

Examination of the three charts in figure 7 also indicate that the 
.maximum negative pressure on the wing occvirs at the leading edge near the 
bovindary between regions 1 and 2. In fact, it can be shown that the max- 
imum occurs at the boimdary. In view of the inequality (66), indentation 
in accordance with the area rule reduces the maximum perturbation veloci- 
ties on the wing. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Oct. 7, 1954 
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Figure 1.- Subsonic pressure coefficient on upper surface of wing of the 

basic combination; mP = 0.5* 
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Figure 2.- Supersonic pressure coefficient on upper surface of wing of 

the basic combination; m3 = 0.5. 
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Figure 3.- Pressure coefficient on 

wing -body 
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upper surface of wing of the indented 
combination. 
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(a) For use with section AA of figures 1 and 2. (h) For use with section BB of figures 1 and 2. 

Figure 5«" Function to be added to give values of pressure coefficient for other Mach numbers. 








